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Pseudogap Phenomena and Superconducting Fluctuations 
in Hubbard Model for High-Tc Cuprates 
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The pseudogap phenomena in High-Tc cuprates are investigated on the basis of the Hubbard 
model which includes only the on-site repulsive interaction U. We consider the pairing scenario 
for the pseudogap. The pseudogap arises from the resonance scattering due to the strong super- 
conducting fluctuations. First, the electronic state and the anti- ferromagnetic spin fluctuations 
are calculated by using the fluctuation exchange (FLEX) approximation. The T-matrix which 
is the propagator of the superconducting fluctuations is calculated by extending the Eliashberg 
equation. The characteristics of the superconducting fluctuations due to the pairing interaction 
arising from the spin fluctuations are represented in the T-matrix. The self-energy due to the 
superconducting fluctuations is calculated by the T-matrix approximation. The pseudogap is 
shown in the single particle properties and the magnetic properties. Moreover, a comprehensive 
explanation of the doping dependence of the pseudogap is obtained. Furthermore, we apply the 
theory to the electron-doped cuprates and obtain the consistent results with the recent experi- 
ments. Finally, the self-consistent calculation for the spin fluctuations, superconducting fluctua- 
tions and the single particle properties are carried out within the FLEX and the self-consistent 
T-matrix approximations. The relation between the superconducting fluctuations and the spin 
fluctuations are clarifled. The calculated superconducting critical temperature Tc is remarkably 
reduced from the results of the mean field (FLEX) calculation. In particular, it is shown that 
the critical temperature decreases with decreasing doping in the under-doped region with large 
U. 

KEYWORDS: High-Tc cuprates; Hubbard model; Pseudogap; Strong coupling superconductivity; 
Superconducting fluctuation; spin fluctuation 



§1. Introduction 

The pseudogap phenomena in under-doped High-Tc 
cuprates have been interested for many years. The pseu- 
dogap is considered to be a key issue for the comprehen- 
sive understanding of the High-Tc superconductivity. 

The pseudogap phenomena mean the suppression of 
the spectral weight above Tc without any long range or- 
der. They are universal phenomena observed in vari- 
ous compounds of High-Tc cuprates in the under-doped 
region. First, the pseudogap was found in the mag- 
netic excitation channel by_the nuclear magnetic res- 
onance (NMR) experiment .Id^ At present, the pseudo- 
gap phenomena have beeiL jobservei in various quan- 
tities which inchjtk NMR,ElB&iB&Q'^eutron scatter- 
ingjs' transport,ElE3' optical spectrum,EJ' electronic spe- 
cific heatjES depsity of statesllS' and the single particle 
spectral weight .ES The experimental results are reviewed 
in ref. 15. 

There are many theoretical scenarios proposed for the 
pseudogap phenomena. An important one is the resonat- 
ing valence bond (RVB) theory. The RVB theory is an 
approach from the non-Fermi liquid state and attributes 
the pseudogap to the singlet pairing of the spinous. t3 As 
an approach from the Fermi liquid state, the magnetic 
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scenario has been investigated. In this scenario, the 
pseudogap is an anti-ferromagnetic gap formation or its 
precursor. 

In this paper, we derive the pseudogap phenomena 
which arise from the strong superconducting fluctua- 
tions. The strong superconducting fluctuations origi- 
nates from the strong coupling superconductivity and 
the quasi-two dimensionality^^ This theory belongs to 
the pairing scenario in which the pseudogap is a precur- 
sor of the superconductivity. Some aspects of the ob- 
served pseudogap have supported the pairing scenario. 
In particular, the measurements of the single particle 
properties have given the important features of the pseu- 
dogap. The same energy scale and the same momen- 
tum dependence between the pseudogap and the super- 
conducting gap have been shown by the angle-resolved 
photo-emission spectroscopy (ARPES) x3 The tunneling 
spectroscopy have also shown the smooth change^ from 
the pseudogap state to the superconducting state.E5 The 
close relation between the pseudogap state and the su- 
perconducting state clearly indicates that the pseudogap 
is a precursor of the superconductivity. 

The paring scenarios have beer 
tively investigated in recent years 
They are classified into some kinds. The scenario based 
on the phasepSjUctuations has been proposed by li^W 
and KivelsonEj) and calculated by other authorsE3'E3 
The strong phase fluctuations in the ordered state are 
expected in the under-doped region since the London 
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penetration depth A is longEJ^ The small London con- 
stant A — I/AttX^ means the small phase stiffness. 

The importance of the strong coupling supercoaduc- 
tivity has been pointed out by Randeria et ax3 on 
the rUcps of the Nozieres and Schmitt-Rink (NSR) the- 
ory.CJ'o' The NSR theory describes the crossover from 
the conventional BCS superconductivity to the Bose con- 
densation of tightly-bound pre- formed pairs. The con- 
cept of the NSR theory is that the crossover is described 
by only shifting the chemical potential /i. In case of 
the strong attractive interaction, the fermions form the 
pre-formed pairs without phase coherence. In this case, 
the gap opens in the spectrum of the fermions. The low 
energy excitations correspond to the pre-formed bosons. 
The NSR theory has been investigated by many authors 
because the gap in the f errn i npic e ycitatinris is consid- 
ered to be the pseudogapBj'OL§E3'E3^ The existence of 
the pre-formed pairs has been proposed phenomenologi- 
cally by Geshkenbein et alx3 with reference to the sign 
problem of the fluctuational Hall effect .£3 

The strength of the superconducting coupling is in- 
dicated by the ratio /e-p- Here, ef is the effective 
Fermi energy, and is the superconducting critical 

temperature obtained by the mean field theory. Because 
the assumption of the weak coupling /e-p <C 1 is jus- 
tified in usual superconductors, the BCS mean field the- 
ory correctly describes the superconducting phase transi- 
tion. The strong coupling superconductivity means that 
the assumption is violated and therefore the supercon- 
ducting fluctuations play an important role. The effects 
of the fluctuations are furthermore strong in the quasi- 
two dimensional systems. The small ep, high-Tc and the 
quasi-two dimensionality are the characteristics of the 
cuprates. Therefore, it is natural to consider the strong 
superconducting fluctuations in High-Tc cuprates. 

The NSR theory takes into account the correction of 
the superconducting fluctuations on the chemical poten- 
tial ^. However, we have asserted that the NSR theory 
is justified in the low density limit and not in High-Tc 
cuprates which are high density systems.H Therefore, 
the NSR scenario should be ruled out as a candidate of 
the pseudogap phenomena. 

In the high density systems, the effect of the strong su- 
perconducting fluctuations manifests in a different way. 
It has been shown that the pseudogap phenomena are 
derived by the self-energy copception due to the strong 
superconducting fluctuations.ll3'E3 The resonance scat- 
tering by the low energy superconducting fluctuations 
gives rise to the anomalous features of the self-energy 
and leads to the pseudogap. These effects have been 
neglected in the NSR theory. The large weight of the su- 
perconducting fluctuations at the low energy necessarily 
exists when the strong coupling supermnd-uctivity occurs 
in the quasi-two dimensional systems Jlj'tiP We adopt the 
resonance scattering scenario also in this paper. 

We furthermore emphasize that our scenario is clearly 
differentiated from the NSR scenario. The phase transi- 
tion in the NSR theory is regarded as the Bose conden- 
sation of the pre-formed pairs. However, the phase tran- 
sition is the supeijCDnductivity with strong fluctuations 
in our scenario.oE^ The later is the realistic situation 



in High-Tc cuprates. 

The effects of the resonance scattering are especially 
weak in the usual weak coupling superconductors. In 
these cases, the self-energy correction has not been em- 
phasized in the theory of the superconducting fluctua- 
tions, since they are less singular compared with the cor- 
rection on the two-body correlation function, such as the 
AL term and the MT term.cZlEa' However, the supercon- 
ducting fluctuations have strong effects on the electronic 
state in the strong coupling case. 

The resonance scattering scenario properly explains 
the magnetic field dependences of the pseudogap phe- 
nomena^ -Lneaaiired by the high field NMR experi- 
ments.oEj'Eil'L3 The calculation gives a comprehensive 
explanation of the experimental results including their 
doping dependenceo-' Moreover, We have explained the 
pseudogap phenomena in various experiments such as 
ARPES, tunneling spectroscopy, NMR, in- plane and c- 
axis transport and the optical conductivityE^ 

Thus, the pseudogap phenomena are well explained 
on the basis of the resonance scattering scenario. How- 
ever, most of the calculations are based on the model 
with an attractive interaction. Our previous calcula- 
tions are also based on the effective model in which the 
effectively strong d-wave paixing, mteiiaction affects the 
renormalized quasi- particles.ll§BE3'BEf Actually, the 
d-wave pairing interaction arises from the repulsive in- 
teraction between the electrons. 

Here, it is one of the most important issues to describe 
the microscopic theory which starts from the repulsive 
interaction and derives the superconducting fluctuations 
and the pseudogap phenomena. It should be confirmed 
that the superconducting coupling is strong enough to 
lead to the pseudogap phenomena in the under-doped 
region. Such microscopic calculation will naturally re- 
produce the doping dependence oLthe pseudogap phe- 
nomena. As is mentioned before,E3' there are two im- 
portant factors in order to realize the strong coupling 
superconductivity. One is that the effective Fermi en- 
ergy e-p is renormalized by the electron-electron corre- 
lation. The other is that the high critical temperature 
Tc is obtained by the pairing intcracticp, mediated by 
the anti-ferromagnetic spin fluctuations .BEj-* The spin 
fluctuations also result from the strong electron-electron 
correlation. Thus, the both factors are derived from the 
electron-electron correlation and should be described in 
a unifled way. The main purpose of this paper is describ- 
ing the microscopic theory by which the superconducting 
fluctuations and the the pseudogap phenomena are de- 
rived from the electron-electron correlation. 

Some authors have calculated the pseudogap phenom- 
ena bypStpjtiag from the model with an repulsive inter- 
action.EJEJ'Ld^ However, the sufficient result describing 
the pseudogap phenomena has not been obtained with- 
out any phenomenology. 

In this paper, we start from the Hubbard model which 
is an typical model for the strongly correlated electron 
systems. First, we describe the quasi-particles and the 
anti-ferromagnetic spin fluctuations by using the FLEX 
approximation. The characteristic momentum depen- 
dence arising from the anti-ferromagnetic spin fluctua- 
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tion£3'EB' are described by the FLEX approximation. 
The superconducting fluctuations are derived from the 
pairing interaction mediated by the spin fluctuations. 
The derived superconducting fluctuations include the 
characteristic momentum and frequency dependence of 
the spin fluctuations. The self-energy due to the super- 
conducting fluctuations is obtained by the one-loop ap- 
proximation (T-matrix approximation). As a result, the 
comprehensive understanding is obtained for the pseu- 
dogap phenomena on the basis of the resonance scatter- 
ing scenario. The doping dependence which includes the 
electron-doped cuprates is also explained consistently. 

This paper is constructed as follows. In §2, we explain 
the Hubbard Hamiltonian and the FLEX approximation. 
In §3.1, we explain the formalism which describes the 
superconducting fluctuations arising from the spin fluc- 
tuations and the self-energy due to the superconduct- 
ing fluctuations. The calculated results for the single 
particle properties are shown in §3.2. In §3.3, we show 
the characteristics of the superconducting fluctuations 
including their doping dependence. The magnetic prop- 
erties in the pseudogap state are explained in §3.4. The 
consistent results with the NMR and the neutron scat- 
tering experiments are obtained. In §3.5, the theory of 
the superconductivity and the pseudogap phenomena is 
applied to the electron-doped cuprates. The relevant re- 
sults including the particle-hole asymmetry are obtained. 
In §4.1, we clarify the relation between the superconduct- 
ing fluctuations and the spin fluctuations in details. In 
§4.2, we show the results of the self-consistent calculation 
including the spin fluctuations, the superconducting fluc- 
tuations and the single particle properties. The critical 
temperature reduced by the fluctuations are calculated. 
As a result, the appropriate phase diagram is obtained. 
In §5, we summarize the obtained results and gives some 
discussions. 



High-Tc cuprates (see Fig. 13). We define the hole- 
doping concentration S — I — n. The main part of this 
paper is concerned with the hole-doped cuprates S > 0. 
A brief application to the electron-doped cuprates 5 < 
will be carried out in §3.5. The nearly half-filling sys- 
tems, which are interested here, should be regarded as 
high-density systems. Therefore, the shift of the chem- 
ical potential, which is considered in the NSR scenario, 
is not important in this paper. The second term in eq. 
(2.1) expresses the on-site repulsive interaction between 
the electrons. All of the phenomena described in this 
paper result from the electron correlation effects. 

The FLEX approximation has been used in order to 
describe the strongly correlateiLelectron systems near the 
anti-ferromagnetic instability.EJ The FLEX approxima- 
tion is a conserving approximationEEP and describes the 
anti-ferromagnetic spin_fliictuations including the mode 
coupling effects partly.Ej'EJ^ It has been shown that the 
superconducting critical temperature with a relevant or- 
der Tc ^ lOOK is obtained-, bjf |-th^ EliEX-approxima- 
tion for High-Tc cuprates.BBByME§E3) The rele- 
vant results are also obtained for the. pcaajiic supcrcon- 



ductor K-(BEDT-TTF) compoundsE 



LJi_|and for the 

ladder type compound Sri4_a:Caa;Cu2404iEl' Moreover, 
the FLEX approximation is used to explain the anoma- 
lous properties resulting from the anti-ferrpmagnetic spin 
fluctuations, such as the Hall coefficient Jl3' 13 the neu- 
tron respnance peak and the dip-hump structure in 
AKPES E3 

The self-energy given by the FLEX approximation 
SF{k,iu!n) is expressed by the one-loop diagram ex- 
changing the normal vertex Vn{q,ifln) (Fig- 1(a))- 



q.in„ 



iri„). 



(2.3) 



§2. Hubbard Model and FLEX Approximation The normal vertex is expressed as follows. 

First, we explain the Hubbard model and the FLEX 
approximation. The Hubbard model has been used for 
a long time as one of the typical models describing the 
strongly correlated electron systems. The Hamiltonian 
is described as, 



k,s 



k,k\q 
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K(q,if^n) = c/^[-Xs(g,if^n) + ipcM^^n) - xo(g,if^«)]- 

(2.4) 

Here, Xs(g,if^n) and Xc(Q',i^^n) represent the spin and 
charge susceptibility, respectively, 

Xa{q,i^n) 



(2.1) 



The first term is the kinetic term. In this paper, we 
adopt the two-dimensional dispersion relation given 
by the tight-binding model for a square lattice includ- 
ing the nearest- and next-nearest-neighbor hopping t', 
respectively, 

Efc — — 2i(cos kx -I- cos ky) + At' cos k^ cos ky — ^. 

(2.2) 

We fix t = 0.5, t' = 0.25i and the lattice constant a — 1. 
In this case, the band width W — A. The chemical po- 
tential fi is determined according as the filling n. These 
parameters well reproduce the typical Fermi surface of 



Xs(g,if^„) 

Xc(9,if^rO 



1 - Uxoiq,i^n)' 
xo{q,i^,i) 



(2.5) 



l + Uxoiq^i^n)' 
where xo(Qii^^n) is the irreducible susceptibility, 

Xo{q, if^n) = X! i^")G'(fc + q, icj„ + if7„). 

k,iuj„ 

(2.6) 

The first order term in Vn{q,ii^n) with respect to U 
is eliminated because it gives only the trivial Hartree- 
Fock term. In the FLEX approximation, the dressed 
Green function G(fc,iu;„) = (iw„ — £k ~ ^F(fc,iw„))~^ 
is used. The self-energy and the spin susceptibility are 
determined self-consistently. Equations. (2.3-6) are self- 



4 



Youichi Yanase, Kosaku Yam ADA 




Fig. 1. (a) The normal self-energy and (b) the anomalous self- 
energy obtained by the FLEX approximation. 



consistently solved by the numerical calculation. In 
the main part of this paper, we divide the first Bril- 
louin zone into 64 x 64 lattice points for the numer- 
ical calculation. The spin susceptibility given by the 
FLEX approximation Xs(q,iri„) is enhanced near the 
anti- ferromagnetic wave vector Q = (tt, tt). The anti- 
ferromagnetic spin fluctuations described by Xs{Q,^^n) 
play a dominant role in the FLEX approximation. The 
characteristic results-. p£ .the . p^adji-piliJeja-omagnetic 
Fermi liquid theoryBS ySB B BBBEI) are qual- 
itatively reproduced within the FLEX approximation. 

The superconducting critical temperature Tc is deter- 
mined as the temperature below which the linearized 
Dyson-Gor'kov equation has a non-trivial solution (Fig. 
1(b)). The criterion for Tc is described by the Eliashberg 
equation which is the following eigenvalue equation, 

X(j){k, iw„) = - T ^ Va(fe - p, iujn - it^m) 

X |G(p,iw™)p0(p,iw„). (2.7) 

Here, Va(q,irin) is the anomalous vertex for the singlet 
channel. This is given by the FLEX approximation as 
follows, 

(2.8) 

The critical temperature Tc is obtained as the temper- 
ature where the maximum eigenvalue A,nax becomes the 
unity Ajnax = 1- The eigenfunction 0max(p, iwm) corre- 
sponding to the eigenvalue Amax is the wave function of 
the Cooper pairs. In this paper, the symmetry of the 
superconductivity is always the dx2_y2-waMe. 

Here, we show the typical results of the FLEX ap- 
proximation in Figs. 2-4. The results for the analyti- 
cally continuated self-energy Sp{k,uj) are shown in Fig. 
2. In this paper, the analytic continuation is carried 
out by using the Fade approximation. The real part 
KeSp {k,uj) shows the negative slope and the imaginary 
part lmSp{k, lu) has the minimum absolute value at the 
Fermi level. They are the characteristic behaviors of the 
Fermi liquid theory which is violated in the pseudogap 



state. 

Some notable features are present in the self-energy. 
The imaginary part shows the oj-linear dependence near 
the Fermi level contrary to the oj-square dependence in 
the conventional Fermi liquid theory. The behavioiulias 
the same origin as that of the T-linear resistivityEJ-ra 
The imaginary part at w = 0, that is the quasi-particle's 
damping, is large near (tt, 0) ('hot spot') and is small 
near (7r/2,7r/2) ('cold spot'). The above momentum de- 
pendence plays a crucial role in the coherent itL-plcpfi 
transport and the incoherent c-axis transport .BEj' £3* 
The renormalization factor Z^^ = 1 — dKeSp{k,uj)/duj 
has the qualitatively same feature as the quasi-particle's 
damping. A lot of the low energy states near (tt, 0) arise 
from the Van Hove singularity and the large renormal- 
ization factor Zf^^. Because the Fermi surface is trans- 
formed by the electron correlation, the Van Hove sin- 
gularity approaches the Fermi level much more. The 
transformation of the Fermi surface is a characteristic 
feature of the nearly anti-ferromagnetic systemsEEP and 
is confirmed by the FLEX approximation.O' The effec- 
tive Fermi energy for the superconducting fluctuations 
are renormalized by the above features. 

The results for the single particle spectral weight and 
the density of states (DOS) are shown in Fig. 3. In spite 
of the unconventional Fermi liquid behaviors, the typi- 
cal single peak structure appears in the spectral weight. 
Thus, the picture of the quasi-particlcs holds in the 
nearly anti-ferromagnetic Fermi liquid. In particular, 
the pseudogap is not seen. The spectral weight is re- 
markably broad at the 'hot spot' and sharp at the 'cold 
spot' reflecting the momentum dependence of the quasi- 
particle's lifetime (see Fig. 2(b)). The low energy states 
near (tt, 0) give the large DOS near the Fermi level. In 
this paper a small constant damping is added in the self- 
energy when calculating the DOS in order to exclude the 
finite size effect. This treatment has no significant effect 
on the obtained results. 

The obtained results for the superconducting critical 
temperature Tc are shown in Fig. 4. It is seen that 
the critical temperature increases with decreasing the 
doping concentration S, or increasing the repulsive in- 
teraction U. In other words, the critical temperature in- 
creases with the development of the anti-ferromagnetic 
spin fluctuations. The critical temperature tends to be 
saturated in the large U and the small S region. It is a 
common result of the FLEX calculation because the spin 
fluctuations have not only the pairing effects but also the 
de-pairing effects. 

§3. FLEX+T-matrix Approximation 

Hereafter, we consider the superconducting fluctua- 
tions and their effects on the single particle proper- 
ties and the magnetic properties. In our scenario, the 
superconducting fluctuations play a dominant role for 
the pseudogap phenomena. Therefore, we have to de- 
scribe the superconducting fluctuations derived from the 
anomalous vertex T4(<7, if^n) given by the FLEX approx- 
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Fig. 2. The self-energy obtained by the FLEX approximation, 
(a) The real part, (b) The imaginary part. Here, U = 1.6, 
S = 0.095 and T = 0.010. The solid lines and the long-dashed 
lines correspond to (Uti", ^'i") ('hot spot') and (^tt, fj''") {'cold 
spot'), respectively. The value TieEp (k, 0) is positive at the 'hot 
spot' and negative at the 'cold spot'. This means that the Fermi 
surface is transformed to be more appropriate to the nesting. 

imation. Although the calculation has been difficult 
we succeed in the description including the momentum 
and frequency dependence of the anti-ferromagnetic spin 
fluctuations. We explain the method of the calculation 
bellow. The obtained results well explain the pseudogap 
phenomena including their doping dependence and some 
characteristic properties. In this section, the couphng 
constant U is fixed to U = 1.6 unless we specify. 

3. 1 Formalism 

First, we explain the formalism of the FLEX+T- 
matrix calculation which includes the effects of super- 
conducting fluctuations within the lowest order. The 
superconducting fluctuations are generally represented 
by the T-matrix which is the propagator of the super- 
conducting fluctuations. The T-matrix is expressed by 
the ladder diagrams in the particle-particle channel (Fig. 



Fig. 3. (a) The single particle spectral weight and (b) the density 
of states obtained by the FLEX approximation. The parameters 
are the same as those in the Fig. 2. 

5 (a)) and is derived from the following Bcthc-Salpeter 
equation. 

T ( fci, iwn : fca, iw,„ : q, iri„) = K(fci - ^2, i^n - i^m) 
- T ^ Va(fei - fc, iw„ - \uji)G{k, iLOi) 

X G{q - k,inn - iuJi)T{k,iLUi : k2,iuj„i : q,iri„). (3.1) 

Here, the pairing interaction is given by the anomalous 
vertex obtained by the FLEX approximation. Generally, 
it is difficult to solve the above integral equation except 
for the case whece-the separable pairing interaction is 
assumed.taEil'ES'ESlE^' Therefore, we carry out the fol- 
lowing two approximations where the meaningful com- 
ponent as the ci2;2_j^2-wave superconducting fluctuations 
is properly taken out. 

The T-matrix at q — f7„ = is approximately decom- 
posed into the eigenfunction with, its respective eigen- 
value of the Eliashberg equation.c^ 

T (fci,itj„ : k2,iuJrn ■ Q = if^n = 0) 
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Fig. 4. The superconducting critical temperature obtained by the 
FLEX approximation, (a) The doping concentration dependence 
for U = 1.6. (b) The repulsive interaction dependence for S = 
0.09. 
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Fig. 5. (a) The T-matrix. (b) The self-energy due to the super- 
conducting fluctuations. 



E 



(3.2) 



The eigenvalue Aq and the eigenfunction (j)a{k,iu>n) 
are derived from the Ehashberg equation, eq. (2.7). The 
index a denotes each mode included in the T-matrix. 
Here, we take out the component with the maximum 
eigenvalue Amax and the corresponding eigenfunction 
^max which has the d^2_y2--wsLve character. By consid- 
ering that the superconducting transition is determined 
by the condition Amax — 1, we can understand that the 
mode described by Amax and (/)max represents the dx^^y2- 
wave superconducting fluctuations. The function 0max 
is the wave function of the fluctuating Cooper pairs in 
the fluctuating regime. Hereafter, we neglect the other 
modes since they have no significant effect on the su- 
perconducting fluctuations and the electronic state at 
the low energy. Hereafter, we neglect the index 'max' 
for simplicity. The approximation is justified when the 
superconducting fluctuations are strong. By using the 
above approximatinn. the T-matrix is expressed by the 
following equations Ej^ 

T(fci,icj„ : k^,iujm ■ Q^^^n) = 
5(q,iri„)</)(fei,ia;„ : q, iri„)(/)* (fc^, i^^m : qA^n) 



1 - A(q, iiln 



(3.3) 



The function A(q,ir2„) is given by the maximum eigen- 
value and (f>{k,iujn '■ q,iQn) is the corresponding eigen- 
function of the following equation. 



A (q,ifi„)0(fc,iw„ : q,il7„) = -T ^ K(fc -p,i 



: q,iO„)p = 1, and 
given as, 5r(q,if7„) = 
q,iri„)Va(fci — fc2,it^« - 



X G(p, icj„i)G(q - p, if7„ - icj„i)(/)(p, icj„ : q, if7„). (3.4) 

The function 1 — A(q,iil„) describes the disper- 
sion relation of the fluctuating Cooper pairs. In 
the above expression, the wave function is normal- 
ized as, X)fe,i^„ 
the coupling constant is 

iu;m)0(fc2 ; i^m ■ q^ i^n)- 

Here, we perform one more approximation, which is 
valid around q = ri„ = 0. We neglect the q- and r2„ 
dependence of the eigenfunction (f)(k,iLUn ■ Qji^^n)- The 
wave function is fixed to that at q ^ fi„ = 0. This ap- 
proximation is well justified in the following way. The 
main q- and ri„ dependence of the eigenvalue A(q,ir2„) 
is derived from the differential of the Green function 
G{q — p,ifln~iuJm) in eq. (3.4). In other words, the dis- 
persion of the fluctuating Cooper pairs, 1 — A(q,ir2„) is 
insensitive to the change of the eigenfunction. ft should 
be noticed that the wave function (j3(k,iLUn '■ q,iQn) is 
independent of q and 11™ when the separable paring in- 
teraction is considered. H 

We have explicitly solved the eigenvalue equation (eq. 
(3.4)). The results have confirmed that the q- and il„ de- 
pendence of the eigenfunction (/)(fc, iujn '■ q, ^n) is small 
in the region where the dominant contribution to the 
self-energy is yielded. In particular, it is confirmed that 
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the value 1 — A(q,iil„) is found within the error of 10- 
20 par cent by this approximation. Thus, the second 
approximation only slightly underestimates the eigen- 
value X{q,i^n) around q = fin — 0, since the precise 
eigenfunction is determined so as to optimize the eigen- 
value. Needless to say, the precise eigenvalue is obtained 
at q = il„ = 0. In other wards, the TDGL parameter 
(see eq. (3.10)) b is overestimated by the second approx- 
imation whereas the mass term to is obtained precisely. 
The approximation slightly underestimates the effects of 
the superconducting fluctuations. Since the approxima- 
tion is not precise for large q and Qn, the unphysical 
contribution to the self-energy is caused by the region. 
Therefore, we eliminate the q- and ri„ independent term 
g(j){ki,iujn)4'* {k2,iuJm) in the T-matrix. 

By using the above two approximations, the T-matrix 
is expressed as follows, 

T (fci,icj„ : k2,iujm ■ q,i^n) 
__ ff A(g, i»„)0(fci , icjn)^* (fca, i^m) ^^ 
1-A(q,if}„) ' ^^-^^ 

where 

A(q,iri„) = -T ^ ^ (?!)*(fc,ia;„)Va(fc — p, iw„ - iwm) 

)0(p,iw™). (3.6) 
Here, the coupling constant g is defined as, 

9^ <l>*i^i^^n) 

k^ ,iiJn k^ jiuJm 

X K(fci — k2,iujn -iuJm)4>{k2,iuJm), (3.7) 
and the wave function is normalized as, 

J2 |^(fe,ic.„)p = l. (3.8) 

By taking out the freedom of the superconducting fluc- 
tuations, we define the pair susceptibility as t{q, if2„) ~ 
g/{l — X{q,iQn)) which_Kp.^]jLav£i.termed the T-matrix in 
the previous papers. Ej'oBlj'EI 

The self-energy due to the superconducting fluctu- 
ations is given by the T-matrix approximation (Fig. 
5(b)). 

, i^.^ : q^ i^n) 

q,i^n 

X G{q - k,in„ -iLJn). (3.9) 

In this section, we use the Green function determined 
by the FLEX approximation GF{k,iuj„) = {itOn — £fc — 
IJF{k,icjn))~^ in calculating eqs. (3.5)-(3.9). That is to 
say, we calculate the lowest order correction due to the 
superconducting fluctuations on the FLEX approxima- 
tion. We call the calculation FLEX+T-matrix approxi- 
mation. The self-energy is obtained by the summation, 
S{k,iu!n) = SF{k,iLUn)+Ss{k,ii^7i)- We show the results 
of the FLEX-I- T-matrix approximation in the following 
subsections. Since the calculation is carried out with the 
fixed chemical potential /i, the doping concentration d 
decreases with decreasing the temperature T. However, 



the difference is smaller than 0.01 and has no significant 
effect. 

In §4, we carry out the self-consistent calculation in 
which the fully dressed Green function G{k, iujn) = 
(iwn — Efc — ^Fikjiujn) — ^sik,itOn))~^ IS uscd every- 
where. As a result of the self-consistency, the effect of 
the superconducting fluctuations are reduced. However, 
the qualitatively similar conclusions are obtained. 

3.2 Pseudogap in the single particle properties 

In this subsection the pseudogap, which is the main 
subject of this paper, is derived. The calculated results 
well justify the pairing scenario based on the self-epergy 
correction due to the superconducting fluctuations.E3'L3' 
On the basis of the model with a d-wave attractive inter- 
action, it has been shown that the anomalous jjroDerties 
of the self-energy give rise to the pseudogap.ll3'Ej' The 
characteristics of the self-energy is that the real part has 
the positive slope and the imaginary part has the max- 
imum absolute value at the Fermi level (see Fig. 8 in 
ref. 19). They are anomalous compared with the typical 
behaviors in the Fermi liquid theory (see Fig. 2). The 
large imaginary part reduces the single particle spectral 
weight at the Fermi level and gives rise to pseudogap. 
The anomalous features seem to compete with the Fermi 
liquid behaviors obtained by the FLEX calculations. In 
particular, the negative slope of the real part which is 
related to the renormalization factor Z^^ is generally in- 
creased by the strong electron correlation. However, the 
calculated results clarify the important mechanism of the 
pseudogap formation in the strongly correlated electron 
systems. 

We show the obtained results for the analytically con- 
tinuated self-energy S^{k,u!) in Fig. 6. Here, the 
doping concentration corresponds to the under-doped 
cuprates. It should be noticed that the Fermi liquid be- 
haviors are seen when we look with the large energy scale 
Lo ~ 0.5. However, the anomalous behaviors leading to 
the pseudogap are clearly seen in much smaller energy 
scale LO ~ 0.05. The anomalous behaviors vanish around 
(7r/2,7r/2) because of the d^2_y2 wave symmetry of the 
fluctuating Cooper pairs. 

It is the important point that the superconductivity 
and the pseudogap take place in the renormalized quasi- 
particles near the Fermi surface. In other words, the 
pseudogap occurs with much smaller energy scale than 
that of the electron systems. The above results clearly 
show the smaller energy scale and justify the calculation 
based on the effective model in which a d-wave attractive 
interaction acts on the renormalized quasi-particles. 

The results for the spectral weight are shown in Fig. 
7. The pseudogap clearly appears in the single particle 
spectral weight (Fig. 7(a)) and the DOS (Fig. 7(b)). 
The dx2-y2-'wave form of the pseudogap is naturally ob- 
tained by the d2.2_j^2-wave symmetry of the supercon- 
ducting fluctuations (Fig. 7(a)). The inset of Fig. 7(b) 
shows the much smaller energy scale of the pseudogap 
compared with the band width. We wish to note again 
that the pseudogap is derived from the self-energy correc- 
tion due to the superconducting fluctuations which are 



8 



Youichi Yanase, Kosaku Yam ADA 



(a) 



0.4 



3 0.2 



CD 
DC 



0.0 



-0.2 _ 









(0.4571, 0.43;u) . 














V V 

















-2.0 



-1.0 



0.0 

CO 



1.0 



2.0 



0.0 



(b) 



^ -0.2 

3 

i -0.4 



-0.6 





/ 


\ (0.98IC, 0.02it) 




/ 


V (0.45JC, 0.437C) 




/ 
1 

/ / 
/ / 






/ / 
/ / 






/ / 
/ / 






/ / 






/ / 
/ / 
// 




" 




— 









-2.0 



-1.0 



0.0 

CO 



1.0 



2.0 



Fig. 6. The self-energy obtained by the FLEX+T-matrix approx- 
imation. (a)The real part, (b) The imaginary part. Here 

U = 1.6, & = 0.095 and T = 0.010. The solid line and the long- 

/63 



dashed line correspond to ( |^ , 
(cold spot), respectively. 



64 



7r) (hot spot) and ( 



64'^) 



enhanced by the strong counling .superconductivity and 
the quasi-two dimensionahty.llSESEJ) The w-dependence 
of the wave function furthermore advances the pseudo- 
gap formation. Thus, the resonance scattering scenario 
is justified by the microscopic calculation based on the 
Hubbard model. It is microscopically confirmed that 
the superconducting fluctuations are sufficiently strong 
to give rise to the pseudogap. 

The energy scale of the obtained pseudogap is consis- 
tent with that of the superconducting gap. The exper- 
imental results show that the energy scale of the pseu- 
dogap is_n£arly the same as that of the superconduct- 
ing gapMEJ The ratio 2AjTc ^ 12 has been obtained 
by thii FLEX approximation in the optimally-doped re- 
gionO^ Here, Ag is the maximum value of the gap func- 
tion in the ordered state. The larger ratio is expected in 
the under-doped region. Our results for the under-doped 
cuprates show the ratio 2Apg/Tc ~ 20 (Fig. 7(a)). Here 
Apg is the energy scale of the pseudogap. Thus, the pseu- 
dogap obtained by our calculation seems to have the rel- 



evant energy scale compared with the superconducting 
gap. 



(a) 



3 
< 



15.0 



10.0 



5.0 



0.0 



" (0.45lt,0.437l) 1 ] 




(0.83it,0.147i) ] 1 




(0.98it.0.02ji) 1 1 








/ ' ' ' 









-0.2 



-0.1 



0.0 

0) 



0.1 



0.2 



(b) 



0.8 



3^ 
q1 



0.6 



0.4 




— FLEX 

— FLEX+T-matrix 



-0.4 -0.2 0.0 0.2 0.4 

CO 



Fig. 7. (a)The single particle spectral weight obtained by the 
FLEX-|-T-matrix approximation. The solid and long-dashed 
lines correspond to those in Fig. 6, respectively. The dash- 
dotted line corresponds to (^ TT, ^vr). (b) The DOS obtained 
by the FLEX (long-dashed line) and the FLEX-|-T-matrix (solid 
line) approximations. The inset is the same result shown in the 
large energy scale. The other parameters are the same as those 
in Fig. 6. 



We show the detailed results for the DOS in Fig. 8. We 
choose the temperature T = 1.25Tc in all figures of Fig. 
8 in order to fix the distance to the critical point. The 
doping dependence is shown in Fig. 8(a). It is shown 
that the pseudogap becomes weak with increasing the 
hole-doping. The gap structure is filled up and the DOS 
near the Fermi level increases in the optimally-doped re- 
gion. The effects of the superconducting fluctuations al- 
most disappear in the over-doped region S > 0.2. This 
is because both the critical temperature and the renor- 
malization of the quasi-particles are reduced by the hole- 
doping. Since the ratio Tc/e-p decreases, the supercon- 
ducting coupling also decreases. Therefore, the effects of 
the superconducting fluctuations are reduced. Thus, the 
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above calculation well explains the doping dependence of 
the pseudogap phenomena. 
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Fig. 8. (a) The doping dependence of the DOS obtained by the 
FLEX+T-matrix approximation. The solid, long-dashed and 
dash-dotted lines correspond to 5 = 0.095, 0.156 and 0.190, re- 
spectively, (b) The density of state for the various strength of U . 
The thin solid, dashed, long-dashed, dash-dotted and thick solid 
lines correspond to C/ = 1.2, 1.3, 1.4, 1.5 and 1.65, respectively. 
In all figures, the temperature is chosen as T = 1.25Tc. 



Here, the importance of the strong correlation is shown 
in Fig. 8(b). It is shown that the pseudogap becomes 
remarkable by increasing JJ . The reason is the same as 
that in the explanation of the doping dependence (Fig. 
8(a)). Thus, the strong superconducting fluctuations 
and the resultant pseudogap are the characteristics of 
the strongly correlated electron systems. 

3. 3 Superconducting fluctuations 

In this subsection, we discuss the character of the su- 
perconducting fluctuations in order to emphasize the im- 
portance of the strong couclmg superconductivity. 

In our pervious papers Jlj'O the TDGL expansion is 



used for the pair susceptibility in the following way. 

to + bq^ - {ai + m2)n' ^ ' 

where the factors arising from the wave function 
(/)(fe, \LOn) are neglected. In the notation of this pa- 
per, the TDGL expansion corresponds to the expan- 
sion for the eigenvalue function A(q, fi) as 1 — A(<7, 51) — 
to + bq^ - (ai + 102)^- 

The TDGL expansion parameters generally describe 
the character of the superconducting fluctuations. The 
detailed properties of the TDGL parameters are dis- 
cussed in ref. 19. Here, we discuss the TDGL parameters 
on the basis of the calculated results in this paper. 

The parameter to = 1 — A(o, 0) represents the distance 
to the phase transition, and is sufficiently small near T^. 
The parameter 02 expresses the time scale of the fluc- 
tuations. The parameter ai is usually neglected within 
the weak coupling theory because it is higher order than 
02 with respect to the superconducting coupling Tc/e-p. 
However, the parameter jii should not be neglected in 
the strong coupling case.EJ^ The calculated results show 
that the sign of ai is negative in the under-doped re- 
gion. The sign of ai which is fixed by the particle-hole 
asymmetry determines the sign of the Hall conductivity 
arising from the superconducting fluctuations .oLlP The 
Hall anomalyca' has been explained by assuming the ex- 
istence of the electron-like pre-formed pairs.ES However, 
it has been pointed out that the explanation is not nec- 
essarily justified within the theory of the strong coupling 
superconductivity.EJ' 

The parameter b represents the dispersion relation of 
the fluctuations and is related to the superconducting 
coherence length as 6 oc ^q. The parameter b is an 
important factor determining the strength of the fluctu- 
ations. The small b generally means the strong super- 
conducting fluctuations. Therefore, we show the results 
of the TDGL parameter b in Fig. 9. 

The TDGL parameter b is expressed by the Fermi 
liquid description within the weak coupling theory as 
b = ffdPd(O) 32^2 Here, gd is the d-wave compo- 
nent of the residual interaction ZfJ/ii[k — k')Zj^i , namely 
5d = Efc,fc"^d(fe)^fcK(fe - k')Zk.c^l{k'). The effec- 
tive density of state of the quasi-particles Pd(e) is de- 
flned as Pd(£) = Sfc'^(^ ~ -E'fc)l0d(fc)P, which is en- 
hanced by the renormalization. Here, is the quasi- 
particle's energy E}^ = Z]^{ek + Rei7^(fc,0)). The vp 
is the effective Fermi velocity for the d-wave symmetry. 
This is defined as Vp = /p pfepi)(/eF)^dfcF/pd(0), where 
the integration is carried out on the Fermi surface and 
Pkp = \(t>d{kF)\'^ /v{kp). The velocity ^(/cf) is the abso- 
lute value of the quasi-particle's velocity = dE/^/dk 
on the Fermi surface. It should be noticed that the quasi- 
particle's velocity is reduced by the renormalization, es- 
pecially around (tt, O). The renormalization is caused by 
the mass renormalization Zf. and the transformation of 
the Fermi surface. The bare velocity Vf, = de^/dfe is 
small at the 'hot spot' on the transformed Fermi surface. 
Since the effective Fermi velocity vp is mainly determined 
by the region around (tt, 0), the velocity vp is remarkably 
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Fig. 9. (a) The doping dependence of the TDGL parameter b at 
the critical point T = Tc- Here, the coupling constant is fixed 
to f7 = 1.6. The inset shows the same quantity at the fixed 
temperature T = 0.0080 (b) The TDGL parameter b for the 
various strength oi U at T = Tc . Here the doping concentration 
is fixed to 5 = 0.09. The inset shows the same quantity at the 
fixed temperature T = 0.0082 



reduced by the electron correlation. 

It should be noticed that the TDGL parameter b is 
proportional to Vp/T'^, which is proportional to the in- 
verse square of the superconducting coupling Tc/ep. In 
this sense, the effective Fermi energy should be defined as 
Ef oc ■Up- That is to say, the TDGL parameter b decreases 
with increasing the superconducting coupling. We can 
see in Fig. 9 that the TDGL parameter b decreases with 
under-doping and/or with increasing U. It is because 
the critical temperature Tc increases and the renormal- 
ization for the effective Fermi velocity vf becomes re- 
markable. These effects reduce the TDGL parameter 
b in spite of the increasing coupling constant (7dPd(0). 
The effects of the quasi-particle's renormalization is con- 
firmed by showing the results at the fixed temperature, 
(see the inset in Fig. 9(b)) Although the couphng con- 
stant gdPd(O) increases and the temperature T is fixed, 
the parameter b decreases with increasing U. This is 



because the effective Fermi velocity v-p is reduced. 

Thus, the superconductivity becomes strong coupling 
one when the electron correlation is strong. The above 
results microscopically justify our scenario about the 
doping dependence of the pseudogap phenomena. The 
superconducting fluctuations become strong, and the 
pseudogap phenomena take place in the wide temper- 
ature region with decreasing the doping concentration S. 
These doping dependence is consistent with the experi- 
mental results. Moreover, the above doping dependence 
of the parameter b is consistent with the magnntip fip1d 
dependence of the pseudogap phenomena.Oc3'L2(EirE3' 

3.4 Magnetic properties 

In this subsection, we show the results for the mag- 
netic properties which are measured by NMR, neutron 
scattering and so on. The anti-ferromagnetic spin fluc- 
tuations are one of the important characters of High- 
Tc cuprates. The development of the spin fluctuations 
are well described by the FLEX approximation. One 
the other hand, the suppression of the low frequency 
spin fluctuations iave been pointed out by NMR mea- 
surements.BB&Q'Er&EP The phenomena have been called 
'spin gap' at the initial stage. At present, it is known 
that not only the spin channel but also thu^aingle par- 
ticle properties show the gap-like featuresll3'EfP There- 
fore, the 'spin gap' is considered to be a result of the gap 
formation in the single particle properties, namely 'pseu- 
dogap'. In this subsection, the pseudogap phenomena in 
the magnetic properties are explained by considering the 
effects of the superconducting fluctuations on the single 
particle properties. 

In the FLEX-|-T-matrix approximation, the spin sus- 
ceptibility Xs{q,^) is obtained by extending the FLEX 
approximation. 



Xo(g,iaO ^-TJ2 G{k,\Lo,,)G{k + g. 



(3.11) 

iri„), 

(3.12) 



where the dressed Green function G{k, luon) — {^^n ~ 
— SpikjioJu) — X's(fc, iw„))^^ is used. The effects of 
the superconducting fluctuations are included in the self- 
energy, Ss{k,iujn)- 

The NMR spin-lattice relaxation rate l/Ti and spin- 
echo decay rate I/T2G are obtained by the following for- 
mula. 



1/TiT 



1 



Q 



Fi(q)[-Imxf(9,c^) 



(3.13) 



1/1:2 



2G = Et^il (9)i^exf (g, 0)]2 - [J2 F\\ (9)Rexf (g, 0)]^ 

(3.14) 

Here, F^{q) = ^[{Ai + 2B{cosq^ + cosgy)}^ -|- {A2 + 
2B{cosqx + cos qy)}"^] and F||(q) = {A2 + 2B{cosqx + 
cosgj,)}^. The hyperfine coupling constants A1-A2 and 
B are evaluated as Ai = 0MB and A2 = -4B.E3' 



Pseudogap Phenomena and Superconducting Fluctuations in Hubbard Model for High-Tc Cuprates 



11 



The calculated results for the NMR l/TiT, \IT2G and 
the static spin susceptibility are shown in Fig. 10. The 
results show the pseudogap in the NMR 1/TiT (Fig. 
10(a)). In the FLEX calculation, the NMR l/TiT in- 
creases with decreasing the temperature (see the inset 
in Fig. 10(a)). This corresponds to the Curie- Weiss low 
of the spin fluctuations. In the FLEX+T-matrix cal- 
culation, the NMR l/TiT increases with decreasing the 
temperature from high temperature, shows its peak at T* 
and decreases with decreasing the temperature. This de- 
crease above the critical point T =:T>^ijajlJii.^_ivcll- known 
pseudogap in NMR measurements.BBIj'QLl'lj tP This phe- 
nomenon is caused by the superconducting fluctuations. 
Since the DOS is reduced by the superconducting fluc- 
tuations, the low frequency spin fluctuations are sup- 
pressed. Because the NMR \/TiT measures the low fre- 
quency component of the spin fluctuations, \/TiT de- 
creases with approaching the critical point. Thus, the 
pseudogap observed in NMR 1/TiT takes place through 
the pseudogap in the single particle properties. 

The NMR I/T2G also shows the pseudogap (Fig. 
10(b)) with the same onset temperature T* as that in 
l/TiT. The NMR I/T2G decreases with approaching 
the critical point. This is also an effect of the super- 
conducting fluctuations. However, the pseudogap in the 
NMR I/T2G is weak compared with that in the NMR 
l/TiT. This is because the NMR I/T2G measures the 
static spin susceptibilityj-jKhich reflects the total weight 
of the spin fluctuations .E^ES It should be noticed that 
the pseudogap suppresses only the low frequency com- 
ponent of the spin fluctuations. This is a natural result 
because the superconductivity has a smaller energy scale 
than that of the spin fluctuations. Thus, it is no wonder 
that the scaling relation of the^uin fluctuations is vio- 
lated in the pseudogap state.d'cI'LrLl-' The above results 
for the NMR give the qualitatively ^ijopidslfint explana- 
tion for the experimental resultsBaou'EtQl]^ Minutely 
speaking, the different behaviors of l/TaCpha-tp been re- 
ported for different High-Tc compounds.oD'Q'EP There is 
an idea that attributes, ihe difference to the effects of 
the interlayer coupling.QQ^ The results of the single layer 
compcjunds show the decrease of I/T2G in the pseudogap 
statcQ^ Anyway, the relatively weak effect of the pseu- 
dogap on I/T2G than on 1/TiT is observed in common, 
which is consistent with our results. 

The simflar features of the NMR l/TiT and liT2G 
have been observed in the superconducting state.u' In 
particular, the I/T2G remains even in the low tempera- 
ture, although the l/TiT rapidly decreases. These fea- 
tures arc-the characteristics of the d-wave superconduc- 
tivity.HEH' Therefore, the above results for the pseudo- 
gap state are natural because the pseudogap is a precur- 
sor of the c?-wave superconductivity. 

While the many quantities shnw the pseudogap with 
the same onset temperature T* i}3 the uniform spin sus- 
ceptibility x^(o, OLdecreases from the much higher tem- 
perature than T* However, the decrease of the uni- 
form susceptibility becomes more rapid near T*& We 
consider that the rapid decrease is caused by the su- 
perconducting fluctuations. The calculated results con- 
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Fig. 10. The temperature dependence of (a) the NMR 1 /TiT and 

(b) the NMR I/T2G calculated by the FLEX-|-T-matrix approx- 
imation, (c) The temperature dependence of static spin suscep- 
tibility x^CQjO) at g = (0,0) (Open squares) and at q = (tt, tt) 
(Closed circles). Here, the doping concentration is fixed to the 
under-doped region 5 = 0.093 ~ 0.103. The inset in (a), (b) and 

(c) shows the same quantities calculated by the FLEX approxi- 
mation. 
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firm the consideration. Wc show the results for the uni- 
form susceptibility xf'{o, 0) and the staggard susceptibil- 
ity x^CQjO) in Fig. 10(c). The staggard susceptibihty 
shows the pseudogap at the pseudogap onset temper- 
ature T* observed in 1/TiT. On the other hand, the 
uniform susceptibility X?(Oj 0) decreases with decreasing 
the temperature from much higher temperature. The de- 
crease becomes rapid near T*. These results well explain 
results of the NMR measurements-cP 

The slight decrease of the uniform susceptibility is 
shown even in the FLEX approximations, (see the inset 
in Fig. 10(c).) Therefore, the decrease of the uniform 
susceptibility is not necessarily attributed to the super- 
conducting fluctuations. However, the superconducting 
fluctuations significantly affect the uniform susceptibil- 
ity and remarkably reduce the quantity near the critical 
point. 

The frequency dependence of the spin susceptibility 
well describes the character of the pseudogap in the mag- 
netic properties. The results for the dynamical spin sus- 
ceptibility x^(q, fi) at q = Q is shown in Fig. 11. The 
real part is suppressed at the low frequency in the pseu- 
dogap state (T = 0.0082 in Fig. 11(a)). Thus, the mag- 
netic order is suppressed by the superconducting fluctu- 
ations. 

The imaginary part has been measured by thejinelas- 
tic neutron scattering and shows the pseudogapE^ The 
calculated results show that the imaginary part is re- 
markably suppressed at the low frequency in the pseu- 
dogap state (Fig. 11(b)). This is the pseudogap phe- 
nomenon observed by the neutron scattering measure- 
ments. On the other hand, the spin fluctuations develop 
in higher frequency region. In other words, the pseudo- 
gap transfers the spectral weight of the spin fluctuations 
from the low frequency region to the high frequency re- 
gion. Therefore, the total weight is not so reduced by the 
pseudogap. These features are consistent with the exper- 
imental resultsu' and with the above explanation for the 
NMR 1/TiT and 1/T2G- It is notable that the pairing 
interaction arising from the spin fluctuations originates 
in the relatively wide frequency region. Therefore, the 
d-wave pairing interaction is not so reduced by the pseu- 
dogap. The interplay between the superconducting fluc- 
tuations and the spin fluctuations will be discussed in 
§4.1. 

In the remaining part of this subsection, we discuss the 
commensurate and incommensurate structure of the spin 
fluctuations. The incommensurability in the Y-based 
compounds YBa2Cu3 06+i5 has been pointed out by the 
inelastic neutron scattering measurements, and has been 
discussed in connection with the stripe phase H The 
stripe phase in the La-based compounds La2-a;Sra; Cu04 
has been measured by the elastic npiitrpri srattprinpr anH 
has been actively discussed lately.BSHy BEIMN) 

It should be noticed that both commensurate and in- 
commensurate structures are obtained within the FLEX 
approximation (Fig. 12(a) and 12(c)). The incommensu- 
rability ^inc is defined by the peak position of the dynam- 
ical spin susceptibility — (7r±(5inc, ti") and (tt, 7r±5inc)- 
Whether the spin fluctuations are commensurate or in- 
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Fig. 11. The frequency dependence of the staggard spin suscep- 
tibihty x?^(Qif^) calculated by the FLEX+T-matrix approxi- 
mation for under-doped cuprates. (a) The real part, (b) The 
imaginary part. The solid lines, the dash-dotted lines and the 
long-dashed lines correspond to T = 0.0082, T = 0.016 and 
T = 0.05, respectively. 



commensurate is determined by the chosen parameters, 
i', (5 and V . One of the general results in our calculation 
is that the incommensurability b\^c increases with the 
doping concentration b. These results are qualitatively 
consistent with the experimental resultsEj-* The detailed 
agreement with the experimental results has been dis- 
cussed by the FLEX calculation in the normal state.Ed^ 

The other general result is that the superconducting 
fluctuations enhance the incommensurability. Although 
the commensurate peak is obtained by the FLEX ap- 
proximation in the under-doped region (Fig. 12(a)), it 
becomes incommensurate owing to the superconducting 
fluctuations (Fig. 12(b)). The incommensurability Jinc 
increases in the FLEX+T-matrix approximation in the 
optimally-doped region (Fig. 12(d)). The above effect 
also originates in the pseudogap formation around (tt, 0). 
The results well explain the experimental results of the 
inelastic neutron scattering. The neutron scattering ex- 



Pseudogap Phenomena and Superconducting Fluctuations in Hubbard Model for High-Tc Cuprates 



13 



(a) under-doped (FLEX) 




(b) under-doped (FLEX+T-matrix) 




(c) optimally-doped (FLEX) 




(d) optimally-doped (FLEX+T-matrix) 




Fig. 12. The momentum dependence of the dynamical spin sus- 
ceptibility Imx^((?, f2) at O = 0.01. The results for under-doped 
cuprates (<5 = 0.095 and T = 0.010) (a) by the FLEX approx- 
imation and (b) by the FLEX-|-T-matrix approximation. The 



periments have shown that the spin fluctuations change 
from incommensurate to_commensurate with increasing 
the measured frequency.tiP Although the measurements 
have been carried out in the superconducting state in 
many cases, the quahtatively similar behaviors as those 
in the pseudogap state are expected in the supercon- 
ducting state. This is because the pseudogap and the 
superconducting gap have the same d^2_j^2-wave form. 
The effect of the pseudogap or the superconducting gap 
disappears when the measurement is done at the higher 
frequency than the energy scale of the gap. Thus, the 
incommensurate peak in the inelastic neutron scattering 
can be explained without any assumption of the charge 
and/or spin spatial modulation. The calculation includ- 
ing the superconducting order is desired for the detailed 
agreement with the experimental results. 

3.5 Electron- doped cuprates 

In this subsection, we apply the above calcula- 
tion to the electron-doped cuprates. It is well 
known that not only the hole-doped cuprates but also 
the electron-doped cuprates such as Nd2-xCej;Clij04_y 
and Pr2_a;Ce2;Cu04_j, are the superconductorso' The 
electron-hole symmetry is expected within the simple 
Hubbard model including only the nearest-neighbor hop- 
ping t. However, some different properties from the hole- 
doped cuprates have been pointed out for the electron- 
doped cuprates. The anti-ferromagnetic ordered state is 
robust against the electron doping rather than the hole 
doping. The relatively low superconducting critical tep- 
perature Tc is observed in the narrow doping rangeO^ 

It has been suggested for a long time that the 
electron-doped cuprates are an s-wave or the other 
node-less superconductor. Some experiments have sup- 
ported the node-less superconductivity. The exponen- 
tial dependence_nf the magnetic penetration depth has 
been reported.ES The absence_of the zero bias con- 
ductance peak is also reporteclE^_.while it should ex- 
ist in the d-wave superconductor .l3' The possibility of 
some nodeJess superconductivity has been proposed the- 
oretically.Ej^ However, some recent experimental results 
support the da;2_y2-wave superconductivity also in the 
electron-doped cuprates. The power law of the mag- 
netic4ienetration deptlo^ and the zero bias conductance 
peakH are shown. The dj.2_y2-wa.ve form of the super- 
conducting gap is directly measured by the ARPES.lj' 
Moreover, the phase sensitive evidence for the d-wave 
supercctojductivity is also reported by the SQUID micro- 
scope.EHI' 

We consider that the hole- and electron-doped 
cuprates should be understood comprehensively by the 
theory including their respective characters. The anoma- 
lous properties of the Hall coefficient are well explained 
for not only hole- but also electron-doped cuprates on 
the basis |Of the same formalism treating the spin fluctu- 
ationsOO-* Therefore, it is natural to expect the same 
mechanism for the superconductivity arising from the 
spin fluctuations. Actually, the instability of the d-wave 
superconductivity mediated_bv_L[i^^pin fluctuations has 
been reported theoretically.OO 
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We consider that the apphcation to the electron-doped 
cuprates is an important test of the theories on High-Tc 
cuprates. Here, it is shown that our theory properly de- 
scribes the important properties of the electron-doped 
cuprates. The electron-hole asymmetric properties are 
introduced by the next-nearest-neighbor hopping term t' 
in our model. The features of the electron-hole asymme- 
try are naturally explained by our calculations. More- 
over, some theoretical expectations are proposed below. 
The experimental verification is desirable. 




Fig. 13. The Fermi surface of the hole- and the electron-doped 
cuprates. The thick solid line shows the Fermi surface of the hole- 
doped case (5 = 0.10. The long-dashed and the dash-dotted lines 
show the Fermi surface of the electron-doped case S = —0.10. 
The long-dashed and the dash-dotted lines correspond to case 
t' = —0.25t and t' = — 0.35t, respectively. The thin solid line 
shows the magnetic Brillouin zone. It should be mentioned that 
this figure shows the non-interacting Fermi surface. Actually, 
the Fermi surface is transformed by the anti-ferromagnetic spin 
fluctuations .Ej' 



The main difference between the electron- and hole- 
doped cuprates results from the shape of the Fermi sur- 
face. The distance from the Fermi surface to the Van- 
Hove singularity is especially important. The Fermi sur- 
face of the electron-doped cuprates is obtained by lifting 
the chemical potential fi. The typical Fermi surface is 
shown in Fig. 13. The Fermi level is lifted from the 
Van-Hove singularity (tt, 0) where the dispersion is flat. 
Therefore, the DOS is rather small in the electron-doped 
cuprates. The small DOS means that the electron corre- 
lation is effectively weak. On the other hand, the nesting 
around (7r/2,7r/2) is enhanced and the tendency toward 
the anti- ferromagnetic order is robust. In the light of the 
d-p model, the carrier is confined in the Cu-site in the 
electron-doped cases, while it is mainly in the 0-site in 
the hole-doped ones. This fact probably contributes to 
the robustness of the anti- ferromagnetic order. This dif- 
ference may affect the parameter of the Hubbard model 
which is an effective model of the d-p model in the metal- 
lic phase. However, the essential difference is included in 
the properties of the band structure since only the vicin- 
ity of the Fermi surface is important for the low energy 
physics. 



Since we define the doping concentration as 6 — I — n, 
it is negative in the electron-doped case. We should note 
that the numerical calculation is not so correct for the 
electron-doped cases as for the hole-doped ones because 
of some reasons. The main reason is the finite size effects 
which are serious owing to the large velocity, the weak 
correlation and the low temperature. The finite size ef- 
fects on the spin- and superconducting fluctuations are 
also serious by the same reasons. Therefore, we divide 
the first Brillouin zone into 128 x 128 lattice points in 
the calculation for the electron-doped cases. 

For example, the difference from the hole-doped 
cuprates appears in the magnetic properties. We show 
the results of the FLEX approximation for the momen- 
tum dependence of the dynamical spin susceptibility in 
Figs. 14 (a) and (b). In these figures, we use the same 
parameters t, t' and U as those of the hole-doped cases. 
A remarkable feature of the spin correlation is that the 
range of the spin fluctuations is narrow in the momentum 
space. The width of the range and the strength of the 
anti-ferromagnetic correlation are reduced by increasing 
the electron-doping. 

It is notable that the FLEX calculation is difficult 
untill the electron-doping \5\ ~ 0.12 because the anti- 
ferromagnetic correlation is too strong. This fact im- 
plies that the anti-ferromagnetic order is robust in the 
electron-doped cuprates than in the hole-doped ones. 
The narrow range of the spin fluctuations means that 
the mode coupling effects are weak in the electron- 
doped cuprates. This favors the anti- ferromagnetic or- 
der. The other factors probably contribute to the ro- 
bustness of the anti-ferromagnetism. For example, the 
anti-ferromagnetic order is robust when the three dimen- 
sionality (or the interlayer coupling) is strong. Moreover, 
the frustration due to the carrier doping is weak when 
the carrier enters the Cu-site. 

Hereafter, we discuss the superconductivity and the 
pseudogap in the electron-doped cuprates. The super- 
conductivity mediated by the spin fluctuations is also 
derived by using the FLEX approximation. The most 
favorable superconductivity is the dr^2_y2-wBNe. How- 
ever, the superconducting critical temperature is very 
low. The critical temperature Tc higher than 0.002 is 
not obtained in case of the parameter t' = 0.25t. This is 
mainly because the pairing interaction mediated by the 
spin fluctuations is weak in the electron-doped case. The 
narrow range of the spin fluctuations in the momentum 
space weakens the pairing interaction. It should be no- 
ticed that the pairing interaction results from the spin 
fluctuations in the wide region around Q = (tt, tt). In ad- 
dition, the small DOS reduces the critical temperature 
much more. Moreover, the Fermi surface in the electron- 
doped case is disadvantageous to the d2,2_y2-wave super- 
conductivity, compared with the hole-doped case. 

We can obtain a higher critical temperature by 
choosing the hopping parameter t' so as to reproduce the 
Fermi surface of Nd2-a;Ce2;Cu04_y with more accuracy. 
Here, we choose t' = 0.35i and U = 2.0. The commen- 
surate spin fluctuations over a wider range are obtained 
(Figs. 14(c) and (d)). The spin fluctuations give rise to 
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(a) 5=-0.12 




(b) 5=-0.15 




(c) 5=-0.10 




(d) 5=-0.13 




Fig. 14. The momentum dependence of the dynamical spin sus- 
ceptibility lmx^{q,^) at Q = 0.01 calculated by the FLEX ap- 
proximation, (a) S = -0.123 and (b) 5 = -0.150 at t' = -0.25t 
and T = 0.01. (c) S = -0.104 and (d) S = -0.130 at t' = -0.35t 



the c?2,2_y2-wave superconductivity. However, the critical 
temperature is low and the doping range in which the su- 
perconductivity occurs is remarkably narrow compared 
with the hole-doped ones (see the inset of Fig. 15). The 
maximum value of the obtained critical temperature is 

Tc = 0.0045. The doping range is (5 = -0.102 0.107. 

The results for the superconductivity ax£ qualitatively 
consistent with the experimental results. O 
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Fig. 15. The obtained critical temperature Tc in the electron- 
doped case. The inset shows the phase diagram including both 
hole- and electron-doped cases. 



We show the momentum dependence of the wave func- 
tion of the Cooper pairs (f>{k,iujn) in Fig. 16. Fig. 
16(a) clearly shows the d-wave superconductivity with 
the node in the diagonal direction. It is a notable differ- 
ence between the electron- and hole-doped cuprates that 
a wave function shows the sharp momentum dependence 
in the electron-doped case (Fig. 16(a)). This is because 
the quasi-particle and the Fermi surface are more clearly 
defined in the electron-doped cases. In the hole-doped 
cuprates, the many low energy states lie around (tt, 0) 
where the quasi-particles are broad. On the other hand, 
the Fermi surface is removed from (tt, 0) and the electron 
correlation is effectively weak in the electron-doped case. 
Therefore, the order parameter has a large value only in 
the vicinity of the Fermi surface. 

The effectively weak correlation is consistent with the 
T-square resistivity in the electron-doped cuprates.o* 
Our calculation shows the w-square dependence of the 
imaginary self-energy lmSp{k,uj) in the electron-dopjed 
case. This is consistent with the experimental resultsEj' 
and implies the T-square resistivity. The incommensu- 
rate spin fluctuations have been regarded to give rise 
to the above momentum deperuieiice of the wave func- 
tion in the electron-doped case.ll23 However, it is actu- 
ally not the main reason. Actually, the spin fluctuations 
are always commensurate in the electron-doped cuprates 
within our calculations. 

Thus, the comprehensive understanding for the phase 
diagram including the particle-hole asymmetry further- 
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more support the (i^2_j,2-wave superconductivity in the 
electron-doped cuprates. 



(a) electron-doped 




(b) hole-doped 




3 k. 



Fig. 16. The momentum dependence of the wave function 
</i(fc, iuin) at ujn = ttT (a) in the electron-doped cuprates 
(5 = -0.105 and Tc = 0.0040) and (b) in hole-doped cuprates 
(<5 = 0.093 and Tc = 0.0080). 



the superconducting fluctuations because of the large b. 
Therefore, the effects are precisely much weaker. 

Thus, our calculation explains the doping dependence 
of High-Tc cuprates including the particle-hole asymme- 
try. The comprehensive understanding including both 
the hole- and electron-doped systems rather supports the 
pairing scenario for the pseudogap. 

§4. Self-Consistent Calculation 

In this section, we carry out the self-consistent cal- 
culation including the spin fluctuations, superconduct- 
ing fluctuations and the single-particle properties. In 
the self-consistent calculation (We call the SC-FLEX-I-T- 
matrix calculation), eqs. (2.3)-(2.8) and (3.5)-(3.9) are 
solved self-consistently where the fully dressed Green 
function G(A;,ia;„) — (ia;„ ~ £fc ~ ^ik,iLUn))^^ is 
used. By the self-consistent calculation, we calculate 
the critical temperature Tc reduced by the supercon- 
diictincr fluctuations As is described in the previous 
papers,^^!!^ there exists the singularity arising from 
the two-dimensionality. The singularity is actually re- 
moved by the weak three dimensionality which surely 
exists in the real systems. Therefore, we determine 
the critical temperature Tc as the temperature at which 
A(o,0) = 0.98 in order to avoid the unphysical singular- 
ity. This determination corresponds to that the dimen- 
sional crossover from the two-dimension to the three- 
dimension occurs around A(o, 0) = 0.98. The finite 
critical temperature is obtained by this operation. The 
method of the determination makes no significant differ- 
ence on the following results. For example, the critical 
temperature is not so different (about 10% lower) even 
if we choose A(o, 0) = 0.99 as the critical point. In par- 
ticular, the similar doping dependence is obtained. Our 
choice of the parameter A(o, 0) = 0.98 is owing to the 
reliability of the numerical calculation. Hereafter, we 
consider the hole-doped cases. 



From the above results, we expect that the pseudogap 
phenomena are not so remarkable in the electron-doped 
cuprates. Because of the low Tc and the weak electron 
correlation, the superconducting coupling Tc/ep is small. 
Actually, our calculation gives the large TDGL param- 
eter 6 = 30 for the parameter set used in Fig. 16(a). 
Therefore, the superconducting fluctuations are not so 
significant as well as in the standard BCS superconduc- 
tor. 

The calculated self-energy shows only the weak ef- 
fects which are the same order as those in over-doped 
cuprates. For example, the damping —lmS^{k,0) at 
T = 0.005 and 6 = -0.104 is about 1/40 of that for the 
under-doped case in Fig. 6. These results indicate that 
the pseudogap arising from the superconducting fluctu- 
ations is weak in the electron-doped cuprates, even if 
it is seen by some experiments. The results ^are consis- 
tent with the experimental results of ARPESo^ and the 
neutron scatteringlH^ which do not show the pseudogap 
in the electron-doped cuprates. The weak effects give 
rise to the slight pseudogap in our calculation. However, 
the finite size effect seriously overestimates the effects of 



Spin fluctuation and superconducting fluctuation 

First, we clarify the relation between the spin fluctu- 
ations and the superconducting fluctuations. The two 
fluctuations complicatedly couple to each other through 
the single particle properties. Here, we digest the impor- 
tant factors of the relations. 

The superconducting fluctuations result from the spin 
fluctuations which act as the d-wave attractive interac- 
tion. The spin fluctuations give rise to the renormaliza- 
tion of the quasi-particles especially around (tt, 0). The 
renormalization reduces the effective Fermi energy ep for 
the c?-wave superconductivity as well as the critical tem- 
perature Tc. Therefore, the strong superconducting fluc- 
tuations are obtained at the reasonable temperature. 

The effects of the superconducting fluctuations on the 
spin fluctuations have been investigated in §3.4. The 
superconducting fluctuations suppress the low frequency 
spin fluctuations and the anti-ferromagnetic order (see 
Fig. 11). A residual question is whether the feedback 
effects of the pseudogap on the spin fluctuations suppress 
the superconductivity or not. 

In order to answer the question, we calculate the 
feedback effect on the critical temperature Tc. We 
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carry out the Modified FLEX (M-FLEX) calculation 
in which the fully dressed Green function is used only 
in eq. (2.6). In the other equations G^(fc,iu;„) = 
{iujn ~ £k ^ ^Fik,iu>n))~^ is used. Thus, the effects 
of the superconducting fluctuations on the spin fluctu- 
ations are included, however those on the single particle 
properties are not included. The eqs. (2.3)-(2.8) and 
(3.5)-(3.9) are solved self-consistently in the M-FLEX 
approximation. The results of the FLEX, M-FLEX and 
SC-FLEX+T-matrix calculations are shown in Table. L 
Here, we determine the critical temperature by the con- 
dition A(o, 0) — 0.98 for an equity. 





FLEX 


M-FLEX 


SC-FLEX-l-T-matrix | 




0.0084 


0.0098 


0.0031 1 


9</'max 


24.18 


14.23 


20.1529 1 


7h 


0.06970 


0.02747 


0.04277 1 


7c 


0.00995 


0.00971 


0.00298 1 







Table I. The comparison among the FLEX, M-FLEX and SC- 
FLEX-l-T-matrix approximations. The critical temperature Tc, 
the effective pairing interaction g0maxi damping at 'Hot spot' 
7h and that at 'Cold spot'7c at T = Tc are shown. The parame- 
ters are U = 1.6 and 5 = 0.083 ~ 0.096. Here, the self-energy in 
the M-FLEX approximation is that due to the spin fluctuations 



In order to understand the results, it is important that 
the spin fluctuations have not only the paring effect but 
also the de-pairing effect. The former is from the rela- 
tively wide frequency region, and the latter is from the 
low frequency component. The pseudogap remarkably 
suppresses the low frequency component, however the 
total weight is not so reduced (Fig. 11). Thus, the 
pseudogap reduces the de-pairing effect rather than the 
pairing effect. Therefore, the higher critical tempera- 
ture Tc = 0.0098 is obtained by the M-FLEX calculation 
where Tc = 0.0084 is obtained by the FLEX calcula- 
tion. In other words, the feedback effects are advanta- 
geous to the superconductivity. Therefore, the feedback 
effects do not suppress the superconducting fluctuations 
and the pseudogap. Thus, only the properties of the low 
frequency component are not sufficient in order to under- 
stand the relation between the spin and superconducting 
fluctuations. 

Needless to say, the pseudogap in the single particle 
properties reduces the critical temperature. The lower 
critical temperature Tc = 0.0031 is obtained by the SC- 
FLEX-l-T-matrix calculation. The de-paring effect from 
the superconducting fluctuations is rather drastic than 
that from the spin fluctuations. (Therefore, the pseu- 
dogap is easily caused by the superconducting fluctua- 
tions.) 

In order to make the above understanding clear, we 
show the quantities 50maxj 7h ^^^'^ 7c in Table. I. Here, 
the quantity 0niax is the maximum value of the order 
parameter (f){k,iujn) and the quantity ff'/'max represents 
the strength of the pairing interaction. The damp- 



ing at the 'Hot spot' 7h = -ImZ'^(fch,0) at T = Tc 
represents the strength of the de-pairing effect. Here, 
fch = (0.987r, 0.027r). In should be noticed that the ob- 
tained value 7h by the M-FLEX calculation has the min- 
imum value in Table. I, although the temperature is 
highest. This fact shows that the de-pairing effect of the 
spin fluctuations are suppressed by the superconducting 
fluctuations. 

The damping at the 'Cold spot'|^n-^|_^ImZ'^ (fee, 0) 
determines the in-plane transport .BEj'H Here, fee — 
(0.457r, 0.427r). The Table. I shows that the value jc 
is not so reduced in the M-FLEX calculation from the 
FLEX calculation. Thus, the feedback effects through 
the spin fluctuations are small around the 'Cold spot'. 
Since the self-energy from the superconducting fluctua- 
tions Sg'{kc,oj) vanishes at the 'Cold spot', we can un- 
derstand that the fffcjct of the pseudogap on the in-plan£ 
transport is small£3'Lj-' as is observed experimcntally.Lr 
The slope of the T-linear resistivity is slightly reduced 
by the feedback effects. The anisotropy 7h/7c increases 
due to the superconducting fluctuations. This is impor- 
tant for the incoherent c-axis transport in the pseudogap 
state becausftJ-he c-axis transport is determined by the 
'Hot spot'.E§0) 

4.2 Results of the SC-FLEX+T-matrix calculation 

In this subsection, we show the results of the SC- 
FLEX-l-T-matrix approximation. The qualitatively simi- 
lar results to the FLEX+T-matrix approximation are ob- 
tained, although the effects of the superconducting fluc- 
tuations are reduced by the self-consistency. It is notable 
that the calculation treating a stronger electron-electron 
interaction U is possible because the superconducting 
fluctuations suppress the anti-ferromagnetic order. Here, 
we choose the interaction U = 2.4. 

In Fig. 17, we show the self-energy calculated for the 
under-doped case S = 0.073. The competition between 
the pseudogap and the Fermi liquid behavior is shown. 
The small imaginary part and the positive slope of the 
real part are obtained in the narrow region around the 
Fermi level. The anomalous features due to the super- 
conducting fluctuations are shown in larger energy scale. 
These features are qualitatively the same as those ob.- 
tained within the model with an attractive interaction.^ 
It should be noticed that these behaviors are given by 
the self-consistency between the superconducting fluctu- 
ations and the single particle properties. The necessity 
of the behaviors has been explained in ref. 19. 

The obtained spectral weight is shown in Fig. 18. The 
peak is shown near the Fermi level at the momentum 
around (tt, 0). This reflects the Fermi liquid behavior 
of the self-energy in Fig. 17. The typical three peak 
structure shown in the under-doped case 6 = 0.073 (Fig. 
18(a)) is a common feature to the reaults based on the 
model with an attractive interaction. tj' It should be no- 
ticed that the weight around the Fermi level is remark- 
ably reduced by the superconducting fluctuations which 
give the large imaginary part of the self-energy around 
the Fermi level. The spectral weight at the low energy re- 
gion shifts to the high frequency region. The suppression 
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Fig. 17. The self-energy obtained by the SC-FLEX+T-matrix 
approximation, (a) The real part, (b) The imaginary part. Here, 
parameters are U = 2.4, S = 0.07.3 and T = 0.004. Tlie solid, 
dash-dotted and long-dashed lines are correspond to (Ht, ^tt), 
(fl'^' M'^) (11'^' M'^)' respectively. 



of the spectral weight near the Fermi level corresponds 
to the pseudogap. 

The small weight near the Fermi level may not be ob- 
served by ARPES because the weight is so small. The 
ARPES measurements have their resolving power about 
the momentum and the energy. The thermal broadening 
arising from the Fermi distribution function /(w) affects 
the ARPES furthermore. Anyway, the suppression of the 
spectral weight near the Fermi surface is a reliable result 
and consistent with the calculated result. 

As is shown in the inset of Fig. 18(a). the usual single 
peak structure is obtained by neglecting the self-energy 
due to the superconducting fluctuations i7s(fe,iw„), 
namely G{k,\ijJn) = {iun ^ £k ^ ^F{k/iujn))~^ ■ The 
spectral weight is recovered and the three peak struc- 
ture vanishes when the momentum leaves (tt, 0) along 
the Fermi surface. Similarly, the effects of the supercon- 
ducting fluctuations are reduced by the hole-doping (Fig. 
18(b)). The single peak structure is obtained in the over- 
doped case S = 0.165 (the inset in Fig. 18(b)). In other 



words, the pseudogap is reduced by the hole-doping be- 
cause the superconducting fluctuations are suppressed. 
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Fig. 18. The single particle spectral weight obtained by the SC- 
FLEX-|-T-matrix approximation, (a) The under-doped case, 5 = 
0.073 and T = 0.004. Here, Tc = 0.0032. The solid, dash-dotted 
and long-dashed lines are correspond to (Iff, (ff ^''^) 

and respectively. In the inset, the solid line shows 

the same result as in the main figure. The long-dashed lines 
show the result by neglecting the effects of the superconducting 
fluctuations Ss{k, iwn)- (b) The optimally-doped case, 6 = 0.119 
and T = 0.004. Here, Tc = 0.0035. The solid and dash-dotted 
lines correspond to (||7r, ^i") and (Ht, ^'i"), respectively. The 



inset shows the result for the over-doped case, <5 
0.0035 and k = (||7r, ^tt). 



0.165, T : 



The above effect of the superconducting fluctuations 
becomes more clear by showing the DOS in Fig. 19. 
The DOS near the Fermi level is reduced by the super- 
conducting fluctuations and the gap structure appears 
in the under-doped case (solid line). It is confirmed by 
showing the result in which Ss{k, iWn) is neglected (long- 
dashed line) that the pseudogap in the DOS is caused by 
the superconducting fluctuations. The pseudogap is sup- 
pressed by the hole-doping, similarly (dash-dotted line). 
Thus, the pseudogap is properly obtained by the self- 
consistent calculation including their doping dependence. 




Fig. 19. The DOS obtained by the SC-FLEX+T-matrix approx- 
imation. The solid and dash-dotted lines show the under-doped 
case (<5 = 0.073 and T = 0.004) and the over-doped case 
{S = 0.165 and T = 0.0035), respectively. The long-dashed Une 
shows the result for the under-doped case obtained by neglecting 
the self-energy Es{k, iu>„). 



Fig. 20. The phase diagram obtained by the FLEX and the SC- 
FLEX-)-T-matrix approximations. The closed circles show the 
results of the SC-FLEX-|-T-matrix approximation for U = 2.4. 
The closed and open diamonds show the results of the FLEX 
and the SC-FLEX-fT-matrix approximations for U = 1.6. 



At last, we show the obtained phase diagram in Fig. 
20. The superconducting critical temperature Tc sup- 
pressed by the superconducting fluctuations is calculated 
by the self-consistent calculation. It should be noticed 
that the suppression of Tc from the mean field value be- 
comes remarkable with under-doping. This is a natural 
result because the superconducting fluctuations become 
strong with under-doping. In other words, the pseudogap 
develops with under-doping and the reduced DOS gives 
the reduced critical temperature. An important result is 
the following; The critical temperature has the maximum 
value at (5 ~ 0.11 and decreases with under-doping in the 
SC-FLEX-|-T-matrix calculation for U = 2.4, whereas 
Tc goes on increasing in the FLEX calculation. In other 
words, the mean field critical temperature Tc*^^ devel- 
ops with under-doping. However, the decreasing Tc in 
the under-doped region is obtained by considering the 
superconducting fluctuations. 

It is notable to write again that the strength of the su- 
perconducting coupling is indicated by the ratio T^^/ep, 
and not by Tc. Since the effective Fermi energy ep de- 
creases, the superconducting coupling becomes strong 
with under-doping in spite of the decreasing Tc. 

We can see from Fig. 20 that the critical tempera- 
ture for U = 1.6 does not decrease with under-doping 
even in the SC-FLEX-|-T-matrix calculation. Thus, the 
strong renormalization of the quasi-particles due to the 
strong electron correlation plays an important role for 
describing the under-doped cuprates. The other impor- 
tant thing is that we can treat the strong correlation 
by considering the superconducting fluctuations which 
suppresses the anti-ferromagnetic order. If neglecting 
the superconducting fluctuations, the system rapidly ap- 
proaches the magnetic order. Because the magnetic or- 



der is suppressed by the superconducting fluctuations, 
the situation with the strong anti-ferromagnetic spin 
fluctuations persists in the under-doped region. This fact 
also contributes to the phase diagram shown in Fig 20. 

The enhancement of the electron correlation effects 
with under-doping is necessary in order to obtain the 
phase diagram in which the critical temperature Tc de- 
creases with under-doping. This doping dependence is 
probably underestimated in the FLEX approximation. 
The limitation of our calculation mainly arises from this 
fact. However, the doping dependence is partially in- 
cluded even in the FLEX approximation. Therefore, we 
can show the under-doped region in Fig. 20 although the 
decrease of the critical temperature Tc is not sufficient. 
This fact implies that the more sophisticated calculation 
which describes the stronger renormalization will gives 
the more decrease of the critical temperature Tc. Any- 
way, the superconducting fluctuations play an important 
role in describing the under-doped region in the phase di- 
agram. 

§5. Conclusion and Discussion 

In this paper, we have investigated the pseudogap phe- 
nomena starting from the Hubbard model. The super- 
conducting fluctuations have been derived microscopi- 
cally from the electron-electron correlation. We have 
succeeded in deriving the pseudogap phenomena by us- 
ing the microscopic calculation. 

The renormalized quasi-particles and the pairing in- 
teraction via the anti-ferromagnetic spin fluctuations 
arc calculated by the FLEX approximation. The anti- 
ferromagnetic spin fluctuations give rise to the strong 
renormalization around the 'Hot spot'. Therefore, the 
effective Fermi energy s-p for the d-wave symmetry is 
remarkably renormalized. Because of the high critical 
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temperature Tc and the renormalized Fermi energy ep, 
the superconducting couphng Tc/sf becomes strong in 
the under-doped region. As a result, the pseudogap phe- 
nomena occurs owing to the self-energy correction by the 
strong superconducting fluctuations. The results sup- 
port the scenario which is describecLwithiji_th£ model 
with a c?-wave attractive interaction.li3'LllHBEJ^ 

An important progress is that the pseudogap is de- 
scribed as a phenomenon near the Fermi surface. Be- 
cause the obtained pairing interaction affects only the 
quasi-particles near the Fermi surface, it is naturally de- 
scribed that the superconducting fluctuations and the 
pseudogap have the small energy scale compared with 
the electron correlation. These results support the pic- 
ture in which the pairing interaction becomes effectively 
strong for the renormalized quasi-particles. 

Moreover, it should be emphasized that the calcu- 
lation in this paper properly describes the doping de- 
pendence of the pseudogap phenomena. The supercon- 
ducting coupling becomes weak with increasing the hole- 
doping, since the Fermi energy increases and the critical 
temperature decreases. Therefore, the superconducting 
fluctuations become insignificant with hole-doping. The 
obtained doping dependence of the TDGL parameter b 
has confirmed the scenario and is consistent with the 
the piagnetic field dependence of the pseudogap phenom- 
ena.N) 

The application of the calculation to the electron- 
doped case gives a f nn si 1 1 n rl m-t;t, a n H hi g with the 



experimental results.!! 



i23'E23) We have obtained 



the drr 



-wave superconductivity with low critical tem- 



peratures within the narrow doping range. It is shown 
that the superconducting fluctuations are weak and the 
pseudogap is not observed in the electron-doped case. 
The difference from the hole-doped case is derived from 
the particle-hole asymmetry of the band structure. The 
comprehensive understanding including the particle-hole 
asymmetry furthermore supports our scenario for High- 
Tc cuprates. 

The pseudogap in the magnetic properties has been 
investigated and explained properly. We have calcu- 
lated the quantities observed by the NMR and the neu- 
tron scattering. The obtained results have naturally ex- 
plained the characteristics of each quantities in the pseu- 
dogap state which are subtly different from each other. 

We have clarifled the relation between the spin fluc- 
tuations and the superconducting fluctuations which are 
complicatedly connected with each other. It is shown 
that the feedback effects of the superconducting fluctua- 
tions on the spin fluctuations enhance the superconduc- 
tivity. Therefore, the feedback effect does not suppress 
the superconducting fluctuations. 

Of course, the superconducting fluctuations them- 
selves reduce the critical temperature. We have cal- 
culated the reduced critical temperature by the self- 
consistent calculation including the single particle prop- 
erties and the spin- and superconducting fluctuations. 
Qualitatively the similar results for the pseudogap have 
been obtained by the self-consistent calculation. The cal- 
culated critical temperature shows the maximum near 
5 0.11, and decreases with under-doping, although it 



keeps increasing in the FLEX calculation. Thus, we have 
succeeded in describing the under-doped region by con- 
sidering the strong superconducting fluctuations in the 
strongly correlated electron systems. 

Finally, it should be stressed that the calculation 
in this paper starts from the Fermi liquid state and 
properly describes the High-Tc cuprates including the 
under-doped region. The comprehensive understanding 
of High-Tc cuprates has been obtained including their 
doping dependence from the over-doped region to the 
under-doped one. 

The important and remained problem is the descrip- 
tion of the under-doped limit. The calculation in 
this paper can not be applied to the limit because 
the FLEX approximation does not describe the Mott 
transition. The FLEX approximation probably un- 
derestimates the electron correlation and overestimates 
the anti-ferromagnetic correlation. Therefore, the anti- 
ferromagnetic correlation may be too strong in the 
under-doped region in this paper. The momentmn inde- 
pendent component of the electron correlation is prob- 
ably stronger and the anti-ferromagnetic correlation is 
less strong in the actual state. The calculation beyond 
the FLEX approximation is desirable in order to describe 
the under-doped region more precisely. 

However, it is notable that the essence of the pseudo- 
gap phenomena in High-Tc cuprates has been explained 
in this paper on the basis of the strong superconducting 
fluctuations. It is a strong evidence for the pairing sce- 
nario that the sufficiently strong superconducting fluctu- 
ations and the pseudogap phenomena are derived from 
the Hubbard model. 
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